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Abstract 

Electron transport through a double quantum dot system is stud- 
ied with taking into account electron-phonon interaction. The Keldysh 
nonequilibrium Green function formalism is used to compute the current 
and transmission coefficient of the system. The influence of the electron- 
phonon interaction, interdot tunneling, and temperature on the density of 
states and current is analyzed. Results show that although the electron- 
phonon interaction results in the appearance of side peaks in the conduc- 
tance at low temperatures, they are disappeared in high temperatures. 

1 Introduction 

The study of transport through systems fabricated from two quantum dots has 
attracted a lot of attention during last decades because of their molecular-like 
behavior [TJ [2j [3l [4j [5J [6j [Tj [S] . Experimentally, quantum dots (QDs) can be in 
a serial or parallel configuration [HIS]- Discreteness of energy levels, charging 
effects, and interdot tunneling in double quantum dot (DQD) systems result 
in novel and interesting phenomena such as: negative differential conductance 
(NDC) [SJ[TU], zero bias anomaly [TTJ [T^], current rectification 15], Pauli spin 
blockade [13], ratchet effect [14], etc. Furthermore, a DQD system can be used to 
measure the lifetime of a singlet-triplet state[TS]. DQD systems are also promis- 
ing candidates for spintronic applications [TB] and quantum computing [IT] [TH] . 

Coupling between the electronic and vibrational degrees of freedom can sig- 
nificantly affect the performance of the systems fabricated from QDs. The 
electron-phonon interaction (EPI) in nanostructures has been extensively stud- 
ied both experimentally and theoretically P l l ^ l ^ l ^ l ^ l2 ^ l ^ l3rJ i r27 1 
|2"81 [25] , Different models are used to investigate the influence of the EPI 
on the transport characteristics of the system such as: the rate equation ap- 
proach [231 130] j the kinetic equation method [3T] , the Keldysh nonequilibrium 
Green function formalism [25] [26] [27], etc. Most work has been done on the 
single QD systems. However, the EPI in the DQD systems is an important and 
interesting subject needing more consideration. 
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In this work, we study the electron transport through a DQD system using 
the Keldysh nonequilibrium Green function formalism |32) . The issue has been 
recently analyzed by means of a rate equation approach |33j . It is assumed 
that the electrons of each dot interact with the vibrational degrees of freedom. 
The influence of the temperature, interdot tunneling, and EPI strength on the 
spectral function of each dot and transmission coefficient is examined. In the 
next section, the model used to describe the system is presented. Sec. 3 is 
devoted to the numerical results and in the end, some sentences are given as a 
summary. 

2 Model and formalism 

We consider two single level quantum dots attached to the metallic electrodes 
and assume that the electrons of each dot interact with a local vibrational mode. 
The Hamiltonian describing the whole system is given as (K = 1) 

H =Y,£ ak cl k c ak + £i4di+t[dld R + H.C}+ (1) 

ak i—L,R 

^[w^a; + Xi[a\ + ai)m] + ^[V^c^d, + H.C] 

i aki 

where c' k (c a k) creates (destroys) an electron with wave vector k, energy e ak , 
in lead a whereas, d\ (di) is the creation (annihilation) operator in the ith 
dot. The third term describes the interdot tunneling and t denotes the interdot 
tunneling strength, uji is the phonon energy in the dot i, while Xi stands for EPI 
strength in the ith dot. We assume that each dot can have up to one electron. 
With respect to the fact that the energy levels of the QD are controlled using 
a gate voltage and, on the other hand, the on-site Coulomb repulsion is order 
of a few meV, the assumption is reasonable. Assuming weak lead-dot coupling, 
the EPI can be eliminated using nonperturbative canonical transformation [34] , 
H = e s ' He~ s , where S — J2i Xi/oJi[a\ — ai]m. Therefore, Eq.([T]) becomes 

H =Y, £ a kci k c ak + £ l dld l +t[xlx R d[d R + H.C}+ (2) 

ak i—L.R 

^Wiaja; + Y^ckiXiC^di + H.C] 

i aki 

where ii = e% — X?/u)i is the renormalized energy level of the zth dot due to 

_^ir t_ .i 

the polaronic shift and, Xi — e ^ is the phonon shift generator oper- 
ator. Since the dot-lead coupling strength is weaker than the EPI strength, 
i.e., V a ki << Xi, Xi is replaced with its expectation value, i.e., Xi —< Xi >= 
e ~9i{2N ph i +i) ^ w h ere g i — (A^/wi) 2 , and N p h } i denotes the phonon population 
expressed as N p h,i — (exp(((3uii) — 1) _1 with /3 = l/fc^T. In the following, it 
is assumed that both QDs have the same shape and size so that the phonon 
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energy and EPI strength are independent of the QDs indexes and are shown by 
wo, and A, respectively. 

In order to compute the current, the Keldysh nonequilibrium Green function 
formalism and wide band approximation are used so that the current can be 
expressed as [35] 

^ J du{Tr[{f L {w)T L - f R ( u )r R )A{u)] + Tr[(T L r*)(iG<(u>))]} (3) 

where f a (oj)(l + exp((cu — /i Q )/3)) _1 is the Fermi distribution function of the 
lead a, and ^ Q denotes the chemical potential of the lead. A(u>) and G < (ui) are 
the Fourier transformations of the spectral function and lesser Green function, 
respectively. The matrix T a describes the tunneling coupling between the ath 
dot and the leads given as 

v l ( r VaT \ R ( a(3T (3y/aT \ 

\y/aT aT J \/3^r /STo ) [V 

where To is a constant, a describes the difference in the coupling of the electrodes 
to different QDs, and j3 stands for asymmetry in the coupling of the QDs to the 
left and right leads. The spectral function is computed as 

A(u) = i{G > {lo) - G < (u)) = i(G r (uj) - G a {cu)) (5) 

where G r ^ a \u>) is the retarded (advanced) Green function of the whole system. 
To compute G </y> \uj) we follow the procedure introduced in Ref. [27]. Because 
the phonon shift generator operator is replaced with its expectation value, the 
Green functions can be separated to the electron and phonon parts as 

G>(t) = -i< d t (t)d](0) >= <5?.(t)e-*<*> (6a) 

G<(t) = i < 4(0)d,(t) >= <5*(i)e-*<-«> (6b) 

where i(j) = 1,2 stands for the matrix elements, G < ( >s> (t) is the dressed lesser 
(greater) Green function resulting from the electronic part of Eq.([2]), and e - *^ 
comes from the phonon part, < X(t)X'(0) >, given as [34] 



e 



"*W = g[N ph (l - e^) + (N ph + 1)(1 - e^)} (7) 
Using the identity e _ *(*) = ^2 n L n e~ lnuJat , G < ^(uj) can be expressed as 

G>(u) =^T L n G>(uj -nu ) (8a) 



G<(u) =^L n G<{uj + nuo) (8b) 



where L n is a function of temperature, T, EPI strength, and phonon population 
given as [27J 



e ~ B g n 



L n = f- for nhO if T = (9) 



Ln = e -9(2N ph+ l) e ^ In{ 2g^J N ph (N ph + 1)) if T^O 
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where /„ is the nth Bessel function of complex argument. G < (w) can be 
computed using the Keldysh equation and Langreth analytical continuation 
as <2<(w) = <^(w)£:<(w)G (fa;), where £<(w) = i[T L f L (cj)+f R f R (ij)}, and 
r a = r a < X > 2 . Indeed, contribution to E < from electron-phonon interaction 
is disregarded due to weak electron-phonon interaction. In order to compute 
the dressed retarded Green function, the equation of motion technique is used. 
It is straightforward to show that the matrix elements of the Green function 
satisfy the following equation 

[w - ^ - £«]<%(w) = 1 + (f + EfiJG^w) (10a) 
[w - e< - tu]Cr a (u) = (t + £u)GUuj) (10b) 

where i = 2 if i = 1, i = 1 if i = 2, and S = — i/2[r£ + r^]. Indeed, the real part 
of the self energy causing slightly shift in the QD energy levels is ignored (wide 
band approximation). With respect to Eg . (fit)]) , the retarded Green function is 
obtained as 

rrri \ 1 — £2 — £22 t + £12 \ h1 \ 

where 

P(u) = (u- h - En)(w - £2 - E 22 ) - (t + S 12 )(i + £ 21 ) (12) 
The Green function of the whole system can be written as [35] 

G£(i) = -t0(t) < {^(^,4(0)} >= <%.(t) < X{t)X\Q) > (13) 

- e(t)G<(t)[< x\Q)x{t) >-< x(t)xHo) >] 

where after Fourier transformation, it becomes 

G r (uj) = ^L n [G r (uj-nwo)-~G < {uj + nu Jo ) + ^G < (uj-nuj )} (14) 

n 

For simulation purposes, uj q is used as energy unit [37], and it is assumed 
that the energy levels of the QDs are degenerate so that the density of states of 
the QDs are identical shown by A(uj). 



3 Results and discussions 

Fig. 1 shows the density of states as a function of energy for different As and 
ts. In fig. la that t = To, there is a delta-like peak in the bonding state (it — t) 
and a Lorentzian peak in the antibonding state (it +t). For weak As, the below 
relation is obtained for A(u>) at the bonding and antibonding states 

A{ei ± t) = 4r — — ^- (15) 

T (1 - a) 4 + 16(1 ± <Ja)H 2 Y Q 
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For strong As, the bonding and antibonding states are farther from each other 
due to the red shift. Furthermore, satellite peaks are well observed in strong 
EPI coming from the phonon emission by electrons (at negative energies) and 
the phonon emission by holes (at positive energies). Note that the probability of 
the phonon absorption is zero in zero temperature because of N p h = 0. As one 
expects, the spacing between the satellite peaks is equal to the phonon energy. 
The influence of the interdot tunneling strength on A(u>) is analyzed in fig. lb. 
At t/To << 1, the bonding and antibonding states are merged. With increase 
of i, these states are well separated. It is interesting to note that although the 
density of states are independent of the chemical potential of the leads in elastic 
transport, it is dependent on the position of the chemical potential of the leads 
in the presence of the EPI according to Eq.flSJ). However, the dependence does 
not affect significantly on A(to). 

The influence of the temperature on A(oj) is plotted in fig. 2. With increase 
of the temperature, the height of the peaks is increased whereas, their width is 
reduced. Such behavior was previously reported about a single level QD [27] . 
The dressing effect is increased by increase of T, so that T becomes smaller and, 
with respect to the fact that the width of the spectral function is related to the 
broadening, increase of the temperature results in the reduction of the width of 
the peaks. In addition, the probability of the phonon absorption is increased 
by increasing T and, as a result, the satellite peaks in A(u>) result from both 
the phonon emission and the phonon absorption. Unlike the elastic transport, 
the spectral function depends on the Fermi distribution of the leads in inelastic 
transport so that height of the peaks is increased by increase of the temperature 
owing to spreading the Fermi function by increase of T. 

Transmission coefficient, T(u>) — Tr(G a T R G r T L ), is shown in fig. 3 in zero 
temperature. There are two main peaks located in the bonding and antibonding 
states and a antiresonance-like behavior due to the destructive interference be- 
tween different pathways through the system. The effect was before announced 
for a coupled QD system [35]. With increase of the EPI strength, the position 
of the bonding and antibonding states is shifted toward left. Furthermore, the 
secondary peaks are well observed resulting from the existence of the new chan- 
nels for the electron transport because of the electron-phonon coupling. The 
probability of the electron transport through these channels is lesser than the 
probability of the electron transport through the main channels. The effect of 
the dot-lead coupling on T(ui) is studied in fig. 3b. Such analysis was before 
done without considering the EPI [38]. In serial configuration (a = 0), the 
transmission coefficient exhibits two Lorentzian peaks centered at the bonding 
and antibonding states, respectively. In addition, the satellite peaks due to 
the phonon-assisted tunneling are clearly observed. With increase of a, the 
antiresonance-like behavior is also observed because of the destructive interfer- 
ence. 

Fig. 4 shows the current- voltage characteristic of the system for different 
temperatures. In zero temperature, the current exhibits step-like behavior so 
that there is a step when the bonding or antibonding state is located inside the 
bias window. More steps are observed in the presence of the EPI because of 
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the phonon-assisted tunneling. The side peaks are well seen in the conductance 
spectrum. With increase of the temperature, the step-like behavior is vanished. 
Indeed, the electrons gain more energy to tunnel through the system owing to 
the thermal excitation. On the other hand, the effects coming from the EPI 
are not observed in the high temperature so that there is only one peak in the 
conductance. 

4 Summary 

The Keldysh nonequilibrium Green function formalism is used to study the elec- 
tron transport through a double single level quantum dot system in the presence 
of the electron-phonon interaction. The behavior of the system is analyzed in 
both zero and nonzero temperatures using the nonperturbative canonical trans- 
formation. The current voltage characteristic exhibits the step-like behavior in 
low temperatures because of the EPI whereas, such behavior is disappeared at 
high temperatures. The influence of the temperature, EPI strength, and inter- 
dot tunneling strength on the transmission coefficient and density of states is 
also examined. 
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Figure captions 

Figure 1: Density of states as a function of energy in zero temperature. 
Parameters are = 0, \i a = 0, a = 0.4, = 1, and T = 0.2. 

Figure 2: A(u>) versus energy for different temperatures. A = 0.5, and 
t = 0.2. Other parameters are the same as fig. f . 

Figure 3: Transmission coefficient as a function of energy. Parameters are 
the same as fig. f . A = 0.5 in fig. b. 

Figure 4: The current versus voltage for A = 0.5 (solid line) and A = 
(dashed line) at kT = (gray) and kT = 1. Inset shows the conductance. 
Parameters are Si = 3, t = 2r = 0.4. 
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Figure 3: 
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Figure 4: 
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